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Abstract 

We discuss two different methods of obtaining "effective 2x2 
Hamiltonians" of the electromagnetic interaction which include rel- 
ativistic corrections. One is the standard Foldy-Wouthuysen trans- 
formation which we compare with the Hamiltonian obtained from a 
direct reduction of the matrix element of the interaction Hamiltonian 
between positive-energy solutions of the free Dirac equation. The two 
approaches are applied to Compton scattering by a proton for which 
a low-energy theorem exists. It is found that the Foldy-Wouthuysen 
Hamiltonian yields the same result as a covariant calculation. This 
is not true for the direct reduction method which will in general lead 
to incorrect results even after restoring the gauge invariance property 
of the Hamiltonian. Furthermore, it is shown that an identification 
of the Z-diagrams of the usual Dirac representation with the contact 
graphs of the Foldy-Wouthuysen representation is incorrect beyond 
the order of the low-energy theorem. 



1 Introduction 

The use of "effective 2x2 Hamiltonians" for relativistic corrections is one of 
the main ingredients of many calculations in low- and intermediate-energy 
nuclear physics. For the electromagnetic interaction of the nucleon such 
correction terms beyond the standard Pauli-Hamiltonian |J are conven- 
tionally generated through a Foldy-Wouthuysen transformation [0, ||, |]. 
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In some applications the Foldy-Wouthuysen transformation is applied to 
a Hamiltonian containing strong interactions as well |5], [], [Tj . The algorithm 
[0, [3], ||] consists of successive applications of unitary transformations of the 
wave functions. It provides a systematic procedure to block-diagonalize the 
Hamiltonian order by order in 1/M, leading to a decoupling of positive- 
and negative-energy states to any desired order in 1/M. As the unitary 
transformations used are in general time-dependent, the old and the new 
Hamiltonian are not unitarily equivalent || |9] . The new Hamiltonian rather 
contains an additional piece involving the unitary transformation and the 
time derivative of the inverse transformation. Nevertheless, such an approach 
automatically leads to the same result for the S-matrix as the standard co- 
variant calculation |TTJ . 

In a second approach the matrix element of the interaction Hamiltonian is 
evaluated between positive-energy solutions of the free Dirac equation. The 
result is then interpreted as the matrix element of some effective Hamiltonian 
between two-component Pauli spinors. The so constructed Hamiltonian is 
then used in conventional old-fashioned time-ordered perturbation theory. 

In ref. |10| we have studied for a generic interaction how the two methods 
compare and when they are expected to give different results. We have seen 
that it is important to investigate processes which require higher than first- 
order perturbation theory. It is the purpose of this work to apply the different 
methods to a specific process, namely low-energy Compton scattering by a 
proton, in order to make more quantitative statements. For several reasons 
Compton scattering provides an ideal process to test the approaches men- 
tioned above. First of all, when expanded in a power series in the frequency 
to of the photon, the two leading terms of the T-matrix are determined by 
a low-energy theorem [11| The coefficients of the series are expressed 
in terms of on-shell properties of the proton, namely its mass, charge and 
anomalous magnetic moment. The derivation of this theorem is based on 
gauge and relativistic invariance. Secondly, it was shown in ref. (|TTJ that a 



covariant calculation of the s- and u-channel Born terms using the interac- 
tion Hamiltonian for a Dirac particle with an anomalous magnetic moment 
as proposed by Pauli @, [13|] reproduces the low-energy theorem. Finally, we 
can address the important question of gauge invariance which will in general 
impose some restrictions on the construction of effective Hamiltonians. 

Alternatively we could have considered pion photoproduction at threshold 
where two different interactions, namely the strong and the electromagnetic, 
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enter and a low-energy theorem exists as well |L4|, [15], [16], O. However, the 



observations for that process are very similar and all the relevant points can 
be made using Compton scattering. 

Our paper is organized as follows. In the next section we will define 
the interaction Hamiltonians which are used in section 3 to calculate the T- 
matrix in second-order perturbation theory. We will discuss the results to 
order u> 2 /M 3 and comment on gauge invariance. Finally, in section 4 we will 
give a short summary. 



2 Definition of the Interaction Hamiltonians 

In this section we will define the different interaction Hamiltonians which 
we will use to calculate the corresponding T-matrix element for low-energy 
Compton scattering to order 1/M 3 . 



2.1 Relativistic Hamiltonian 

For the relativistic description we will work with the Dirac equation describ- 
ing the interaction of a Dirac proton having an anomalous magnetic moment 
(k = 1.79) with an external electromagnetic field, 

d^f(x) 

t^^ = (H + H I (t)Mx), (1) 

where H = a ■ p + f3M is the free Dirac Hamiltonian and Hj(t) is given by 
1,0 

Hl (t) = (3(eA^-^a^Fn, (2) 

with F^ u = d u A >1 — d^A u . A calculation involving the model Hamiltonian 
of eq. (0) obviously cannot take all strong interaction effects into account, 



such as, e. g. , off-shell effects [18], [19], |(J or transitions to excited states. 
However, for our purposes this is not required, as we only want to use the 
second-order covariant calculation as a reference, which defines the "correct" 
relativistic result within the framework of the model Hamiltonian. 



3 



2.2 The Foldy— Wouthuysen Hamiltonian 

The Foldy-Wouthuysen transformation fj, ^, |J provides a method which 
order by order in 1/M decouples the upper from the lower components. We 
apply the procedure to the Hamiltonian of eq. (|I]) to obtain a representation 
of the 4x4 Hamiltonian which is block-diagonal to order 1/M 3 |[L0|| . We 
define the upper left-hand block of the 4x4 Foldy-Wouthuysen Hamiltonian 
as the effective 2x2 Hamiltonian, H e ^~ FW , to be used in our calculation 
of Compton scattering by the proton in second-order perturbation theory. 
We keep interaction terms up to second order in the coupling constant e and 
obtain, 

i—^— = (H JJ +H 1 JJ (t) + H 2 1J (*))#(», (3) 

with 
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(p 2 A 2 + A 2 p 2 + Ap ■ AA-p+ (V • A) 2 + B 2 



8M 3 

-2ip- AV ■ A + 2zV • A A ■ p + 2p ■ A a ■ B + 2a ■ B A ■ p 



e 2 K 



SM^ v 

+A- V x B-A-e), (4) 



where £ = -V$ - A and B = V x A and £> refers to dO/dt. Note that 
in eq. (|3|) denotes a two-component wave function. 
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2.3 Direct Pauli Reduction Hamiltonian 



In many applications a different method is used to obtain an effective 2x2 
Hamiltonian. Given a general relativistic interaction, it consists of evaluating 
the matrix element of the interaction operator between free positive-energy 
solutions of the Dirac equation and reducing it to two-component form. The 
result may then be expanded in a power series in 1/M and is interpreted as 
an effective operator to be used with a free Hamiltonian Hq^ ~ p . Applying 
this procedure to the interaction Hamiltonian of eq. (|2]) and keeping only 
terms to order 1/M 3 yields the following Schrodinger equation 

with Hq^~ p the same as Hq^~ fw in eq. (|3j) and 

H e / f - p (t) = e<£-—(p-A + A-p)-—(l + K)a-B 
1 y ' 2M KP PJ 2M K J 

[a -p,[a- p, $]] - i-prjA? -p } a-E] 



]M 2 " 11 AM 2 

e r ^2 r -> -> -> X\ i i e r ^2 



-{p\ {a ■ p, a ■ A}} + — — [p\ [a ■ p, a ■ A]] 



gM 3 r> 16M 3 ' 

+ I6^ {a! - P "'' {c? - P "'' (? ^ }} - (6) 

The 2x2 effective interaction Hamiltonian Hj^~ p of eq. (^) may be in- 
terpreted as the upper left-hand block of the transformed operator Hf = 
TqHjTq, where T denotes the time-independent free Foldy-Wouthuysen 
transformation (see e. g. eq. (8) of ref. ||10|| ). We emphasize that Hf is not 
block-diagonal, i. e. , it will connect positive- and negative-energy eigen- 
states of the free Foldy-Wouthuysen Hamiltonian. 

The above procedure only produces a linear interaction term and no con- 
tact interactions. However, it is well-known that a gauge-invariant coupling 
through a minimal substitution in Hq^~ p generates linear and quadratic 
(and higher-order) terms in A as required by gauge invariance. As such 
terms are absent in the Hamiltonian of eq. (|5[) it will not exhibit the correct 
transformation behavior under a gauge transformation. Of course, we can 
cure this problem by introducing additional interaction terms by hand, such 
that the resulting Hamiltonian transforms properly. However, we have to 
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keep in mind that one of our main purposes in considering Compton scatter- 
ing is a critical examination of the use of effective Hamiltonians. In fact, for 
interactions other than the electromagnetic, there is in general no guiding 
principle which tells us whether or not the simple matrix element reduction 
is a reasonable procedure to generate an effective Hamiltonian. For that 
reason we will nevertheless use the Hamiltonian of eq. @ (being aware that 
this will lead to unacceptable results) and will postpone the question of gauge 
invariance until we discuss our results. 



3 Calculation of the T— Matrix Element 

In this section we will calculate the T-matrix element for low-energy Comp- 
ton scattering in second-order perturbation theory using the different in- 
teraction Hamiltonians of section 2, eq. (0), (f|) and ([]), and compare the 
results. 

The kinematical variables and polarization vectors for Compton scatter- 
ing are defined in figure 1. The calculations are performed in the laboratory 
frame, where the energies u and u' of the initial and final photon are related 
by 

= w (l + ^(1 - cos(0))) _1 = u (l-£(l- oo8(6l)) + [^]\ . (7) 

In eq. (0) 9 denotes the angle between the momenta of the initial and the final 
photon in the laboratory system and [uj 2 /M 2 \ means "of order of u 2 /M 2 " . 
Furthermore, we make use of the radiation gauge, $ = 0, V • A = 0. We 
express the result as an operator which still has to be evaluated between 
two-component Pauli spinors of the proton. Any calculation of the T-matrix 
element^ to order 1 / M 3 may be written in the following form [ Yl , 21 ] 



= e ■ e'Ai + id ■ e' x eA 2 + k' x e' ■ k x eA 3 + ia ■ (k' x e') x (k x e)A A 
+ik ■ e' a ■ k X eA 5 + ik! ■ ea ■ k' x e'A 6 

+ik ■ e' a ■ k' x eA 7 + ik' ■ ea ■ k x e'A 8 , (8) 



1 The T-matrix element tfi is related to the invariant amplitude M of ref. j^j through 
tfi = -i(M/^Ef~)M. 
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where k, e(k', e') refer to the direction and the polarization of the initial (final) 
photon, respectively^]. In the following we will calculate the amplitudes Ai 
as predicted to second-order perturbation theory using the different Hamil- 
tonians. 



3.1 Covariant Calculation 

In the covariant calculation the invariant amplitude to second-order pertur- 
bation theory using the Hamiltonian of eq. (0) reads 



M 



-ie 2 u( Pf ) ( £ - — k' A ( A + — k A ufa) 

+ (e^ e',k~ -k% (9) 



where the expression for the u-channel diagram may be obtained from the 
s-channel diagram through crossing symmetry as indicated in eq. @. 
In eq. (^) we have introduced — i8 instead of the standard —it as the small 
imaginary mass required by the Feynman-Stiickelberg boundary condition 
in order to avoid confusion with the polarization vector e M . 

For our discussion it turns out to be useful to split the Feynman propa- 
gator of the nucleon, Sf(p), into its positive- and negative-frequency contri- 
bution, 

S F (p) = - 



b — m + id 
= S { F + \p) + S ( i\p) 

2E(p° -E + i6) 2E(p° + E-iSy 1 ' 

with E = y/p' 2 + M' 2 . Such a separation allows one to identify Z-diagrams 
and ordinary diagrams of old-fashioned time-ordered perturbation theory 
(see e. g. ref. [^3] and figure 2). It should, however, be noted that, although 



the individual diagrams of covariant perturbation theory are Lorentz scalars, 
this is not separately true for their positive- and negative-frequency con- 
tributions f23|. In other words, these contributions are frame-dependent 



2 It is possible to construct a generalized Hamiltonian which parameterizes the T— 
matrix including all strong interaction effects to order uj 2 /M 3 |2l], 22 1. Such an effective 
Hamiltonian involves the Compton polarizabilities a and /3. 
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quantities. Furthermore, as was pointed out in ref. [53], the notion of Z 
diagrams depends on the interaction as well as the representation used. Here, 
we will refer to Z-diagrams (ordinary diagrams) as those resulting from the 
negative- (positive-) frequency part of the propagator of eq. ([[(]) when using 
the Hamiltonian of eq. (@). 

After inserting eq. ([[(]) into eq. @ we expand the resulting expression in 
1/M and bring it into the form of eq. (§). The predictions of this straight- 
forward but tedious calculation for the amplitudes A; separated into their 
forward and backward contributions are tabulated in table 1 and 2. 



3.2 Foldy— Wouthuysen Calculation 

The calculation involving the effective Foldy-Wouthuysen Hamiltonian of 
eq. (f|) is performed in second-order time-ordered perturbation theory and 
the contribution to the S-matrix of second order in the coupling constant 
reads 



T eff-FW\frFW > 
I Oi 



S fi = -2m5(E f + J - Ei - u) j < <$> F f w \H e 2 

^ fi>FW\TTeff-FW\j b FW foFW \ rrej } -FW \foFW 

,3 ^ ^0/ \ n lem I^Op ^0p [-"labs \^0i > 



+ E / d 3 P 

spins 



spins 



Ef + uj'- E^+i5 

^ fr.FW\TrefJ-FW\fcFW foFW I zref f-FW i^FW ^ 

3^< $ 0/ \ H labs \®0p >< ®0p \fhem I^Oi > 

Et — uj — Eg + i5 



(11) 



In arriving at eq. (fTT| ) we have assumed a harmonic time dependence for the 
interaction, i. e. , exp(— iut) and exp(zu/t) for the initial and final photon ap- 
pearing in the interaction Hamiltonian of eq. (^), respectively. Furthermore, 
H\lm FW an d HlH~ FW refer to the emission and absorption of a photon, re- 
spectively, and it is understood that the part of H2^~ FW which at the same 
time creates and annihilates a photon is used. The Hamiltonian H2^~ FW 
gives rise to a contact term, whereas Hl^~ FW contributes "quadratically" 
in second-order perturbation theory (see figure 3). In Coulomb gauge, the 
leading-order term of }j^f~ FW i s Q f order 1/M. As we are interested in the 
T-matrix element to order 1/M 3 , we need the linear interaction Hamiltonian 



8 



to order 1/M 2 onlyQ. Establishing the graphical rules corresponding to the 
interaction Hamiltonians of eq. @ is straightforward, but the expressions 
are lengthy. We only report the predictions for the amplitudes A4 in table 1 
and 2, separated according to their origin as contact or non-contact terms. 



3.3 Calculation Involving the Effective Hamiltonian 
jjyf-P an( j Q aU g e Invariance 

Finally, a second-order calculation involving the effective Hamiltonian of 
eq. (H), Hj^~ p , gives rise to 

S fi = -2iri5(E f + J - Ei - u) x 

^ foFW I tje ff-PlfoFW ~, ^ (hFWlTTeff-PlfoFW . 
^ ^0/ \ n Iem I^Op >< - ^0p \ n Iabs \_^0i > 

E } + uj' - Eg + id 

3 ^ °/ \ n Iabs I^Op ^0p \ n lem I^CH > I 




spins 



Ef — lo — Ef+ id 



(12) 



As we have seen in ref. [1(| such a calculation amounts to neglecting the 
negative-frequency part of the covariant calculation. Thus we already know 
the results for the A{, as we have split the covariant calculation into its 
forward and backward part. 

Approximations, such as the neglect of the negative-energy contributions, 
will in general be in conflict with gauge invariance. Thus we have to address 
the gauge invariance property of the equation of motionF] 

t^^ = H(<f,AMx). (13) 

Gauge invariance means that 

#'(x) = exp(-iex{x))^{x) (14) 

3 If one wants to check gauge invariance to order 1 /M 3 one has to retain terms to order 
1/M 3 in Hl^~ FW , since the leading-order term in jj^ff-FW j s proportional to $. 

4 Although we only write down the potentials $ and A as the argument of the Hamil- 
tonian, it is understood that it may as well depend on derivatives of the potentials. 
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is a solution of 

^o_p — h(® + x,A — VxM*), (15) 



provided ty(x) is a solution of eq. ( |T3| ) [p5[| . From this we may derive the 
following constraint for the form of the Hamiltonian 

H ($ + X , A - Vx) - exp(-iex)#($, A) exp(iex) - ex = 0. (16) 

The effective Foldy-Wouthuysen Hamiltonian of eq. (§) satisfies the gauge 
invariance property to order e 2 whereas the effective Hamiltonian of eq. (|5|) is 
not gauge-invariant even to first order in e. Suppose we denoted the effective 
Hamiltonian of eq. (|5]) by Ha and we wanted to add a Hamiltonian Hb such 
that the sum satisfied the condition of eq. (0). Inserting Ha + Hb into 
eq. ( JT5| ) it is easily seen that the condition for H B reads 

exp(-ie X )H B ( y <5>, A) exp(iex) -H B ($ + x,A- V*) = 

H A ($ + x,A-Vx)~ exp(-iex)H A (<S>, A) exp(iex) - e X - (17) 



From eq. (|17D it is seen that H B cannot be determined uniquely as one may 
always add another term He which separately satisfies 

exp(-iex)Hc(<5>,A)exp(iex)-H c (<5> + x,A-Vx) = 0. (18) 

Let us for convenience consider eq. (^) to order 1/M 2 only, as this is 
sufficient for the point we want to make. Inserting Ha = Hq + H e /^ P 
into eq. flTT|) we obtain the following constraint for Hb 



exp(-iex)H B ($, A) exp(iex) - H B {§ + X,A- Vx) 



e « 2 



*'P,[*-ft *]] + 7777(2^ • V X - (Vx)^ 



8M 2L ^ JJ 2M 



s u2 [X, [* " P, [<? " P, *]]] + j^[x, [B • P, <? • E\\ + [e% (19) 

From the right-hand side of eq. ( |19"D we can now see that H^~ P + Hj 
does not satisfy the constraint of eq. ([Tj]) even to leading order in e. The 
introduction of e. g. 

Hb ^ A) = m& ■ * [ " ■ £ $]] + w + i£ [c? • ^ 5 • ^ ] (20) 
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makes the sum Ha + Hb transform according to eq. (ITJj) to order e 2 and 
1/M 2 . The second term of eq. ( pUP "completes" the minimal substitution in 
the kinetic-energy term of eq. ([5]) to the order we are considering here. The 
other two terms are clearly not unique and a comparison with the Foldy- 
Wouthuysen Hamiltonian to the same order, namely 1/M 2 , shows that the 
two Hamiltonians are in fact different. 



3.4 Results and Discussion 

In table 1 we show the results for the amplitudes to the order of the 
low-energy theorem, namely [1/M] and [uj/M 2 ]. We find that the negative- 
frequency contribution of the covariant calculation (Z-diagrams) yields ex- 
actly the same result as the contact interaction of the Foldy-Wouthuysen 
transformation. Similarly, we have as well a correspondence between the 
positive-frequency part (ordinary diagrams) on the one hand and the non- 
contact terms of the FW calculation on the other hand. 

The naive - but within the framework of nonrelativistic second-order 
perturbation theory consistent - calculation involving the effective Hamilto- 
nian of eq. (^|) is identical with the positive-frequency contribution of the 
covariant calculation and leads to completely unacceptable results. Even the 
leading-order term, the so-called Thomson limit, is not reproduced. Inspec- 
tion of table 1 tells us that the neglect of the negative-frequency contribution 
is responsible for such behavior. Such an approximation will spoil the gauge 
invariance property of the model. However, it is well known that the low- 
energy theorem for the Thomson limit is based on gauge invariance. For 
that reason we constructed an additional interaction term to make the to- 
tal Hamiltonian gauge-invariant. This construction naturally involved the 
e 2 ^4 2 /2M term which then generated the correct Thomson limit. However, 
the other terms are by no means unique. As we can see from table 1 even 
including the additional Hamiltonian Hb of eq. ( PU| ) does not reproduce the 
correct low-energy theorem prediction beyond the leading-order term, al- 
though the Hamiltonian used is gauge-invariant to order e 2 . It was pointed 



out by Low [12| that the determination of the terms beyond the the Thomson 
limit require more input than gauge invariance alone. Our example clearly 
is in agreement with this observation. In other words, enforcing gauge in- 
variance is required to obtain the Thomson limit, but it does not predict the 
next-to-leading-order term of the low-energy theorem. 
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The Compton scattering by a proton firstly demonstrates how mislead- 
ing it may be to use effective Hamiltonians obtained from a two-component 
reduction of a relativistic matrix element. The failure to reproduce the con- 
tributions of the Z-diagrams originates from the fact that the 2x2 effective 
Hamiltonian Hj is generated from a 4 x 4 Hamiltonian H[ which is 
not block-diagonal. Use of H e /^~ P then amounts to a truncation of the 
Hilbert space, namely the neglect of negative-energy states in this approach, 
which will in general lead to unreliable results. In contrast, the 4x4 Foldy- 
Wouthuysen Hamiltonian H FW is block-diagonal and hence using the 2x2 
Hamiltonian H e ^~ FW does not make a difference, as H FW does not connect 
positive- and negative-energy states in the FW representation. 

Furthermore, we have seen that gauge invariance alone is in general not 
sufficient to overcome this shortcoming, i. e. , to determine more than the 
leading-order term. For a similar discussion of the second point within the 
framework of a quantum field theoretical model see ref. |26| . 

The terms of order u 2 /M 3 are beyond the prediction of the low-energy 
theorem and hence are model-dependent. The results of the covariant cal- 
culation and the Foldy-Wouthuysen calculation for these terms are listed in 
table 2. Once again we see that the total result of both calculations is the 



same at this order in agreement with our observations in ref. [10|. However, 



from table 1 one might have concluded that the correspondence between Z- 
diagrams and contact interactions and equally between ordinary diagrams 
and non-contact terms holds true for higher orders in 1/M as well. It can 
now be seen that for some of the amplitudes this is no longer true. In or- 
der to be specific, we find that the Z-diagram contribution to the amplitudes 
Ai, A§ and Aq differs from the contact terms. 

As in table 1 the result of the effective Hamiltonian is given by the 
positive-frequency part of the covariant calculation. We have not made the 
effort of forcing the effective Hamiltonian to be gauge-invariant to order e 2 
and 1/M 3 since the result will not be unique anyway. 



4 Summary and Conclusions 

In this work we considered Compton scattering by the proton within three 
different approaches. Firstly, we performed a covariant calculation to second- 
order perturbation theory. In order to identify Z-diagrams and ordinary 
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diagrams in the covariant approach, we split the Feynman propagator into its 
positive- and negative-frequency contribution. We expanded the so obtained 
T-matrix element to order uj 2 /M 3 and used the result as a reference for the 
nonrelativistic calculations with effective Hamiltonians. 

The calculation involving the Foldy-Wouthuysen Hamiltonian correctly 
reproduced the reference result of the covariant approach. This was expected 
as the 4x4 Foldy-Wouthuysen Hamiltonian is block-diagonal and thus does 
not couple between positive- and negative-energy states. The neglect of the 
negative-energy states in going from the full S-matrix to the effective one 
has no consequence and hence the full (and as a reference correct) S-matrix 
is still reproduced. However, we saw that an identification of the Z-diagrams 
with the contact interaction of the Foldy-Wouthuysen Hamiltonian is not 
correct beyond the order of the low-energy theorem. 

Finally, the effective Hamiltonian H e /^~ p obtained from a direct reduc- 
tion of the relativistic interaction Hamiltonian between positive-energy so- 
lutions of the free Dirac equation led to incorrect results. The reason is 
that Hj^~ P is derived as the upper left-hand block from a 4 x 4 interac- 
tion Hamiltonian which is not block-diagonal. In this case, the neglect of 
the negative-frequency solutions makes a considerable difference. Even the 
leading-order term is not reproduced. This may be cured in part by enforcing 
the gauge invariance property of the corresponding Hamiltonian. However, 
as we saw, such a procedure is not unique, and thus there will always be 
some degree of arbitrariness involved, when enforcing gauge invariance. 
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Table 1: Compton Scattering amplitudes 

Contribution to the amplitudes Ai to the order of the low-energy theorem 
in units of e 2 . The amplitudes A%, Aj and As are zero at this order. The 
second and third column refer to the negative- and positive-frequency con- 
tribution (Z-diagrams and ordinary diagrams) of the covariant calculation at 
second-order perturbation theory. The positive-frequency part of the covari- 
ant calculation is identical with the result of the effective Hamiltonian result- 
ing from a direct Pauli reduction. The fourth column contains the additional 
part due to Hb in order to make the effective Hamiltonian gauge-invariant. 
The fifth and sixth column refer to the contact and non-contact terms of the 
Foldy-Wouthuysen calculation in second-order perturbation theory. 
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Table2: Compton Scattering amplitudes 

Contributions to the amplitudes Ai beyond the low-energy theorem to order 
uj 2 /M 3 in units of e 2 . The amplitudes A2 and A4 are zero at this order. Note 
that we have not constructed a Hamiltonian H B . 
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